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Ft ■ Abstract. We discuss some recent progress in constructing analytic ap- 

proximations to the galaxy clustering. We show that successful models can 

^Sj ' be constructed for the clustering of both dark matter and dark matter 

haloes. Our understanding of galaxy clustering and galaxy biasing can be 
greatly enhanced by these models. 
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1. Introduction 
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Q"^ . The large scale structure of the Universe is believed to have developed from 

small perturbations (usually assumed to be Gaussian) of the matter density 
Q^' field by gravitational instabilities. Under these assumptions the clustering 

pattern and velocity field observed today are determined by the initial 
;-H ' conditions via the perturbation power spectrum [-P(A;)] and the cosmological 

^ . parameters such as fio, the cosmic density parameter. It is therefore possible 

to derive constraints on model parameters from the observed density and 
velocity distributions of galaxies. 
^ . There are two fundamental problems to be addressed: First, since the 

H I galaxy distribution may be biased relative to the mass density field, we need 

to understand such bias before making meaningful comparisons between 
models and observations. Second, even if the observed galaxy distribution 
traces the matter distribution, we still need to understand how the observed 
distribution is related to a cosmogonical model. This is by no means trivial, 
because the clustering pattern and velocity field observed today are non- 
linear. N-body simulations are usually invoked to find a solution to these 
problems. However, such simulations are limited both in resolution and in 
dynamical range, and can be difficult to interpret. Our understanding of 
the underlying physics can be greatly enhanced by simple physical models 
and the analytic approximations they provide. 



In this article we summarize some of our recent progress in connection 
to the problems mentioned above. We show that successful semianalytic 
models can be constructed for the clustering of both dark matter and dark 
matter haloes. Such models can enhance significantly our understanding of 
both the nonlinear evolution of galaxy clustering and galaxy biasing. 

2. Cosmogonies 

The models present below are for CDM-like cosmogonies. The cosmology 
is described by the cosmological matter density (^o), the cosmological con- 
stant (Ao) and the Hubble constant [Hq = 100/ikms^^ Mpc^ ). The initial 
power spectrum is 

P{k)(xkT^{k), (1) 



ln(l + 2.34g) p^ , , ,,_ ,,,,^,2 , ,. .. ^3 , ,« v-, ^4 
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^(^) = 90/ 1 + 3-89g + {IQ.lqY + {5A6qy + (6.71g)^ 

(2) 
where q = A;/(r/iMpc"^), T = Qoh (Bardeen et al. 1986). The RMS mass 
fluctuation in top-hat windows with radius R, cr{R), is defined by 

a\R) = / ^^\k)W\kR), (3) 

JO k 

where W{x) is the Fourier transform of the top-hat window function, and 

A2(fc) = (l/27r2)fc3p(fc) (4) 

is the power variance. We normalize P{k) by specifying erg = cj(8 /i~^Mpc). 

3. Low-Order Statistics of Dark Matter Distribution 

3.1. MASS CORRELATION FUNCTION 

The evolved two-point correlation function ^(r) is related to the evolved 
power variance A|;(fc) by 

^, , f°° dk ^2 ,,, sinfer ,^, 

Thus, in order to get ^(r) we need an expression for AE{k). Following 
the original argument of Hamilton et al. (1991), Jain, Mo & White (1995), 
Padmanabhan et al. (1996) and Peacock &; Dodds (1996, hereafter PD) have 
obtained fitting formulae which relate A^' to P{k) for a given cosmological 
model. The latest version of such a fitting formula is given in PD. The 
solid curve in Fig.l shows the prediction of such a formula for ^(r) in the 
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SCDM: n=l. r=0.5, CTb=0.63 
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Figure 1. Predicted two-point correlation function of mass (solid curve) compared with 
simulation results (symbols). The dashed curve is given by the linear power spectrum. 
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a) open: r=0.a; o-g=1.0_ 




b) open: r=0.a, 0-3=0.5 . 



Figure 2. The redshift-evolution parameter of the mass correlation function as a function 
of redshift, for open cosmologies with fio = 0.1, 0.2, 0.3, 0.5 and 1 (from bottom up). 



standard cold dark matter (SCDM) model with f^o = 1, Aq = 0, /i = 0.5 
and cjg = 0.62 (see Mo, Jing & Borner 1997, hereafter MJB, for details). 
Comparing it with the N-body results, we see clearly that the fitting formula 
works well. This is true for various other cosmogonic models studied. 



3.2. REDSHIFT EVOLUTION OF MASS CORRELATION FUNCTION 

The redshift evolution of the two-point correlation function is usually parametrized 
by the form 

e(r,z)=e(r,0)(l+z)-(3+^-), (6) 

where ^(r, z) is the amplitude of the two-point correlation function at phys- 
ical radius r at redshift z, €p is a parameter to describe the time evolution. 
If ^(r) oc r~'^ , then ep = 7 — 1 for the linear growth in an Einstein-de Sitter 
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SCDM: n=l. r=0.5. crg=0.62 
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Figure 3. Mean pairwise peculiar velocity of dark matter particles. Curve is the model 
prediction; symbols are the simulation results. 



universe, ep = 3 — 7 for clustering patterns fixed in comoving space, and 
ep = for stable clustering (i.e. clustering patterns fixed in proper coordi- 
nates). Given the model in §3.1, it is straightforward to obtain ep. Fig. 2 
shows Ep as a function of z for open cosmologies with various VLq. The value 
of ep at z is obtained by fitting £^{r = l/i~-^Mpc, z') in the redshift interval, 
z' = ^ z. The evolution is more rapid (i.e. ep larger) for universes with 
larger VLq, since linear structures grow faster in a high-f^o universe. For 
high f^O) the evolution is faster at z ~ 1, because of nonlinear evolution. At 
late time when clustering on r ~ l/i~^Mpc becomes stable, the evolution 
becomes slower. 



3.3. MEAN PAIRWISE PECULIAR VELOCITIES 

From the pair conservation equation (Peebles 1980, §71), the ensemble (pair 
weighted) average of the pairwise peculiar velocity {vi2{r)) = ([v(x) — v(x-|- 
r)] • r) can be written as 



{vi2{r)) ^ 1 1 dC{y,a) 

H{a)r 3[l+^(y,a)] 3 In a 



(7) 



where r is the proper, and y the comoving, separation between the pairs; H 
is the Hubble's constant at expansion factor a; ^(y, a) = (3/?/^) j^ y'^dy^{y, a) 
Thus, to obtain {vi2{r)), we need to work out dAE{k,a)/da. This can be 
done directly from the fitting formula of AE{k,a) (see MJB). Fig. 3 shows 
the comparison between the model prediction and the simulation results 



for the SCDM model. The agreement between the two is remarkably good, 
and this is true for many other cosmogonic models studied. 

3.4. COSMIC ENERGY EQUATION 

The (density weighted) mean square peculiar velocity of mass particles 
(vf) is related to the two-point correlation function by the cosmic energy 
equation: 

—a^{vl) = AirGpa^—-^^, 8 

da oma 

where p is the mean density of the universe, and 
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h{a) = j ydyay,a)=J — ^],' ' . (9) 

Integrating eqn. (8) once, we have 

1 



{vj) = ^n{a)H\a)a^l2{a) 



a/2 (a) 



l2{a)da 




(10) 



In the linear case, /2(a) oc D^{a), where D{a) = ag{a) [g{a) is the linear 
growth factor; cq = 1] and 



(vj) = ^n{a)H\a)a^l2{a) 



'-J^)iy'^'^'{ 



(11) 



We found that eqn. (11) is a good approximation (to an error of < 10%) to 
eqn. (10) for all realistic cases. Thus, for a given cosmogonic model, we can 
easily obtain {v\). For various cosmogonies, the model predictions fit the 
simulation results to an accuracy better than 10 percent (see MJB). 

3.5. PAIRWISE PECULIAR VELOCITY DISPERSION 

The relative velocity dispersion of particle pairs of separation r is defined as 
([v(x) — v(x + r)]2)"'^". In Fig. 4 we show (by symbols) the dispersion of the 
pairwise peculiar velocities projected along the separations of particle pairs 
[(■^12 (''))"'^] i^ ^^^ N-body simulations. The main features of {vi2{f))^''^ 
are (a) monotonic rise at small r; (b) saturation at large r; (c) a maximum 
at medium r. As shown in MJB, these features can all be explained by 
physical arguments. 

Based on the N-body results we make the following ansatz for (t'i2(^))' 

{vl2{rt'^ = ^''-''HrMr/rc), (12) 
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SCDM: n=l, r = 0.5, 0-3 = 0. 62 
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Figure 4- Pairwise peculiar velocity dispersion of dark matter particles. Curve is the 
model prediction; symbols are the simulation results. 



where (/'(x) is a universal function and Vc is a nonlinear scale. We choose re 
to be the virial radius of M* haloes, r*. For a given power spectrum, the 
linear radius of M* haloes, Tq , is given hy a{rQ) = 1. The relation between 
r* and Tq in different cosniological models can be obtained analytically, as 
discussed in detail in MJB. We approximate the functional form of (f){x) 
by: 

.( , ^oojl + Bx-^) + Ax-" 

^^''' l + y[5(a)]0-35[^(a)]0-2x-("+«)' ^ ^ 

where 4>^ = ^Jlivl)^/"^ /{Hrc^^-% A, 5, F, a > /3 > 0, and k > are 

constant. {vi2{r)) is forced to have the large separation asymptotic value 

for uncorrelated pairs. For r ^ 0, (wi2(^)) oc x'' so that it increases with 
r as a power law. The solid curve in Fig. 4 shows the prediction of our 
fitting formula with 



A = 58.67; B 



-0.3770; T/ = 4.434; 



a = 2.25; (3 = 1.90; k = 0.15. 
The fit to the simulation data is reasonably good. 



(14) 
(15) 



4. Spatial Clustering of Dark Matter Haloes 

So far we have discussed the clustering properties of dark matter. To com- 
pare model predictions with the observed galaxy distribution, we also need 
to understand how the galaxy distribution is related to the dark matter 



distribution. The bias of the galaxy distribution relative to the mass distri- 
bution can be obtained once we know how galaxies form in the mass density 
field. However, some progress can still be made before the details of galaxy 
formation is understood. In the standard scenario of galaxy formation, a 
gravitationally dominant dissipationless component of dark matter is as- 
sumed to aggregate into dark matter clumps (dark matter haloes), galaxies 
then form by the cooling and condensation of gas within these dark haloes 
(e.g. Kauffmann 1997). It is therefore important to approach the problem 
of galaxy biasing by first understanding how dark matter haloes are dis- 
tributed relative to the mass. In the following we show that simple analytic 
models for such relations can be constructed (see Mo & White 1996, MW). 

4.1. BIAS RELATION 

We define dark matter haloes as virialized clumps of dark matter. In the 
spherical collapse model, a dark matter halo is characterized by its mass, M, 
and the redshift, z, when it is assembled. To describe the relation between 
halo and mass distributions we define a bias relation, 

6h{R) = b{R,6,M,z)6{R), (16) 

where S{R) = [p{R) — 'p]/'p is the overdensity of matter in a sphere of 
radius R, 5h{R) is the same overdensity for dark matter haloes. Suppose 
the conditional density of dark matter haloes within a sphere of radius R is 
nii{M, z\S,R). (This is the number density of dark matter haloes within a 
sphere of radius R, given that the mean mass overdensity within this sphere 
is 5.) The bias relation can then be written as 

A (M .U m nh{M,z\5,R) 

5,iM,z\5,R)= _^^^^^^ -1, (17) 

where nh(M, z) is the mean number density of haloes in the universe, given 
by the Press-Schechter formalism (Press & Schechter 1974, PS): 



_ [^ P ^z da 

n,(M,.)dM = -^---^— exp 
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2^2 ( 



dM, (18) 



where M = {Att / 2,)-pr^ , b-, = 6cD{ao)/D{a), 5c ~ 1.686. The conditional 
density, nh{M, z\6,R), can be obtained by an extension of the PS formalism 
(e.g. Bower 1991): 

/,^ ir T^x,,^ /2"p i^z — 5o)(7(r) da(r) 




Figure 5. Bias relation Sh{5) in a scale free model with n = —0.5. Results are shown 
for spheres with radii R/L — 0.02, 0.05 and 0.13 (L: the side of the simulation box). The 
results for R/L = 0.05 and 0.13 are shifted by 1 and 2 decades along the horizontal axis. 
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Figure 6. Predicted halo-halo two-point correlation functions (solid curves) compared 
to those from simulations (symbols). 



where ^o is the hnear mass overdensity of spherical regions with Lagrangian 
radius Rq. Assuming spherical collapse, we have i?o = [1 + 5{R)Y'^R, and 
5q is determined by S{R) and R through the spherical collapse model (see 
MW for details). With these the bias relation is fixed. Fig. 5 shows the bias 
relation for haloes in a scale free model with P{k) oc k^^'^ . It is clear that 
our simple model works reasonably well. 

4.2. TWO-POINT CORRELATION FUNCTION OF HALOES 
When (5 <C 1 and r <C i?0) the bias relation can be written as 



5h{R;M,z)=b{M,z)5{R)- b{M,z) = l + ^ 

Oz 
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(20) 



Thus the bias factor b depends only on M and z. Under the assumption of 
linear bias, 

Ur)=h\M,z)i{r), (21) 
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a) open: r=0.a; crg=1.0 



Figure 7. The redshift-evolution parameter of the halo-halo correlation function for 
haloes selected according to the first selection rule (see text). Results are shown for open 
cosmologies with fio = 0.1, 0.2, 0.3, 0.5 and 1 (from dot-dashed to solid curves). 
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Figure 8. The same as Fig. 7 for haloes selected according to the second selection rule. 



where ^(r) is the two-point correlation function of mass. Fig. 6 shows the 
two-point correlation functions for haloes in a scale free model with P{k) oc 
^-0.5^ The simple model works reasonably well. 

4.3. REDSHIFT EVOLUTION OF HALO CORRELATION FUNCTION 
As in §3.2, we parametrize the redshift evolution of ^/^ by 

a(r-,^)=a(r,0)(l + z)-(3+^'^). (22) 

Unlike the mass correlation function ^^ the evolution of ^h depends on how 
dark matter haloes are selected. We consider two different selections. In the 
first, haloes are selected at the same time when the correlation function is 
calculated. This case is relevant for galaxies, if galaxies merge as fast as 
their dark haloes. In the second, haloes are selected at an earliear epoch 
than when their correlation function is calculated. This case is relevant 
for galaxies, if they remain distinct after their haloes merge. Fig. 7 shows 
e/j as a function of z in open cosmologies with two choices of cxg. Here 
haloes are selected according to the first selection. As before, the value of 
e/j at z is obtained by fitting ^h{r = l/i~^Mpc, z') in the redshift interval 



10 

z' = — > z. In all cases, the halo correlation function evolves less rapidly 
than the mass correlation function (see Fig. 2), because haloes with fixed 
circular velocities are more biased at higher redshifts. The evolution is also 
less rapid for a lower erg, because of the higher degree of bias involved. 
Fig. 8 shows the same results for haloes selected according to the second 
selection. In this case, the evolution of $,h is faster than that of ^, because 
the haloes (with Vc = 50 — lOOkms"^) are antibiased relative to the mass 
at the time of selection {z = 2). For massive haloes which are {positively) 
biased relative to the mass, the evolution of their correlation function is 
slower than that of the mass. The evolution is also slower for a lower cg, 
because of the increased degree of bias (or decreased degree of antibias). The 
results in this subsection suggest that one needs to be very cautious when 
interpreting the redshift evolution of galaxy correlation function. Without 
knowing in detail what population the observed galaxies are, it is difficult 
to infer the time evolution of the mass correlation from the correlation 
functions of these galaxies. 

5. Discussion 

The models presented above show how the low order statistics of the den- 
sity and velocity distributions are determined by cosmogonies. Thus they 
can be used to construct statistical measures of the density and peculiar 
velocity fields to constrain cosmogonic models by observations. In partic- 
ular, the models can help in the reconstruction of cosmogonic parameters 
from measurements in redshift space. There are also other applications. 
MJB discuss the dependence of the small scale pairwise peculiar velocity 
dispersion on the presence (or absence) of rich clusters of galaxies. Mo, Jing 
& White (1996) use the model to study the correlation of galaxy clusters. 
Mo, Jing & White (1997) extend the model to high-order correlations of 
dark matter haloes and density peaks. 
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